Macroscopic loop formation in circular DNA denaturation 
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The statistical mechanics of DNA denaturation under fixed linking number is qualitatively differ- 
ent from that of the unconstrained DNA. Quantitatively different melting scenarios are reached from 
two alternative assumptions, namely, that the denatured loops are formed in expense of 1) overtwist, 
2) supercoils. Recent work has shown that the supercoiling mechanism results in a BEC-like picture 
where a macroscopic loop appears at Tc and grows steadily with temperature, while the nature of 
the denatured phase for the overtwisting case has not been studied. By extending an earlier result, 
we show here that a macroscopic loop appears in the overtwisting scenario as well. We calculate its 
size as a function of temperature and show that the fraction of the total sum of microscopic loops 
decreases above Tc, with a cusp at the critical point. 

PACS numbers; 87.15.Zg, 36.20.Ey 



The melting or denaturation of DNA refers to the sepa- 
ration of the two complementary DNA strands, a process 
which can be induced either thermally py] or mechani- 
cally dl. The associated phase transition is well under- 
stood by means of theoretical models, e.g., by Poland and 
Scheraga [5|[, Peyrard and Bishop and their various 
extensions [7|-l9| . Thermal melting of DNA forms the ba- 
sis of the PGR technique, while statistical and dynamical 
properties of denatured loops may turn out to be relevant 
for understanding DNA-protein interactions 13, HI and 



gene expression initiation [12|, [13 1 . 



Interestingly, it is a common practice to use plasmid 
(circular) DNAs during PGR, since most bacteria come 
with circular DNA as a means of protection against 
degradation. The resulting entanglement of the two 
strands due to the natural twist of the DNA molecule 
imposes an obvious obstacle for the denaturation process 
(as well as replication, protein synthesis, etc.), which is 
overcome in nature by means of special DNA manipulat- 
ing proteins Nevertheless, the thermal behavior of 
a circular DNA chain in absence of such helper proteins 
proves to be a nontrivial problem and has been addressed 
recently [l5l-[l7|. 

The circular geometry entails the presence of a new 
topological invariant in the system: the number of times 
two chains of the DNA wind around each other, namely, 
the linking number (LK). The thermodynamics of the 
system should therefore be investigated within the cor- 
responding, restricted phase space. This framework is 
also relevant to single-molecule experiments on DNA in 
which the chain ends are rotationally constrained [l8| . 
We here discuss the implications of LK conservation on 
thermal melting characteristics within the framework of 
the Poland-Scheraga (PS) model. 

A consequence of fixing the linking number is that, 
the denatured loops form in expense of (right-handed) 



torsional stress on surrounding DNA duplex segments. 
As with any elastic ribbon with finite bending and twist 
modulus, dsDNA responds to torsion by "supercoiling" 
(bending the backbone as in coiling telephone cords) 
and/or by "overtwisting" (modifying the stacking angle). 
Thermodynamics of a fixed-LK DNA chain whose bound 
segments are unbendable but have finite twist rigidity 
was investigated by Rudnick & Bruinsma fl5'|. The al- 
ternative extension of the PS model considering the pos- 
sibility of supercoil formation, but not overtwisting, has 
also been discussed recently 17, 3]. It has been shown 
that the transition is a kind of Bose Einstein condensa- 
tion (BEG) where the macroscopic loop formed above the 
melting temperature plays the role of the condensate. 

We show below that a BEG-like transition takes place 
in the "overtwisting" scenario, too, which is the main 
contribution of this article. We next compare this phe- 
nomenon with a similar observation we made earlier for 
the supercoiling response and with the denaturation of 
DNA with free ends. We conclude that, the birth of a 
nontrivial macroscopic loop at the melting point is the 
defining characteristics of the thermal denaturation of 
DNA under fixed linking number, irrespective of how the 
molecule responds to torsional stress. 

Ref. [20] gives a detailed account of the melting transi- 
tion in the PS model. The partition function of the model 
can be expressed in closed form and its singular behav- 
ior, with proper treatment of the loop entropy, yields a 
first-order melting transition. A similar analysis is given 



in Ref. [171 . Il9l | for the case of nonzero supercoil density 
under the constraint that the total length of the dena- 
tured loops is proportional to that of supercoils (mimick- 
ing LK conservation). This system, unlike the PS model, 
displays a continuous melting transition, accompanied by 
a loop "condensate" that appears at Tc and grows grad- 
ually with temperature. 
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Overtwisting, i.e., increasing the stacking angle be- 
tween the successive base pairs, is the alternative to su- 
percoiling by which a partially denatured circular DNA 
chain can accommodate the resulting torsional stress on 
duplex regions. An extension of the PS model with over- 
twist has been investigated earlier in [l^, which we ex- 
tend below. In particular, we show here that the melting 
transition in this scenario, too, is accompanied by the 
formation a macroscopic loop. 

To this end, following [l5| , let us consider an arbitrary 
configuration of a DNA chain of a total length L com- 
posed of denatured segments with total length Li and 
bound segments with = L — Li. We assume that the 
linking number expelled by the loops is uniformly dis- 
tributed along the chain and results in a uniform increase 
59 (X Li/Lf, in the stacking angle measured between the 
two ends of a unit DNA segment. Associated overtwist 
energy per unit length is assumed harmonic and the total 
internal energy can be written as 



n 



^2 



(1) 



where k > is a measure of overtwist stiffness in units of 
energy and e;, < is the binding energy per unit length. 
The canonical partition function for the DNA chain can 
then be expressed as [ist 



i2 Lb+lLi + 1 
I z, z, 



where the grand sum Q''^^ — [l/ujZb—l—A^c{szi)]^^ fol- 
lows from the usual PS model with the Boltzmann weight 
uj for a unit bound segment and fugacities {zb,zi} per 
unit length of bound and denatured DNA, respectively. 
The contour integral has a simple pole for Zb which by 
Cauchy formula yields 



Z''{Lb,Li) = 
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F{zi,mi) 



= -(1 - mi) log 



uj[l + A<i>c{szi)] 

^2 



+mi log zi + (3k, 
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where mi^b — Li^b/L and to/ -|- nib = 1- In the thermo- 
dynamic limit, the partition function can be evaluated 
using the saddle-point condition dF = 0. Therefore, 
F{zi,mi) serves as a free energy functional for the DNA 
chain. Minimization yields a continuous phase transi- 
tion for c > 2 governed by the singularity of the polylog 
function at = 1 [l5| . It is straightforward to show 



the critical temperature Tc shifts linearly with the ra- 
tio of the overtwist penalty k and the binding energy 
£6 = -ksTloguj: 
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Here T^^ is the critical temperature of the original PS 
model and to^ is the critical loop fraction which we find 
to be independent of the twist stiffness (therefore equal 
to the corresponding value in the PS model) : 



1 + A(Cc + Cc-l) 



(5) 



Here Cc = ^c{^) is the Riemann zeta function. 

Our goal is to investigate the existence of a macro- 
scopic loop for T > Tc m the above picture. Let us 
assume that such a loop exists with size Lq = fn^L and 
calculate too. As for a Bose gas below the condensa- 
tion temperature, the estimated amount of denatured 
DNA, when calculated as a sum over microscopic loops, 
is now short of the actual value by Lq . Therefore we set 
Li = LJ'""'" + Lo and substitute in EqU] 



Z'^{Li,Lb) — > Z'^{Li — Lo, Lb 



T ^ 



which, following the same steps, now yields an additive 
macroscopic loop correction to Eq|31 



F{zi,mi,mo) = F{zi,mi) - molog{szi) 



(6) 



Note that to/ is the total loop density, including micro- 
scopic and macroscopic contributions. 

For T < Tc {zi < 1/s), the free energy is minimized 
wrt Too at the extremal value too — 0. Therefore no 
macroscopic loop exists below Tc- The total fraction of 
denatured bases (all due to microscopic loops) can be 
calculated by setting dmiF — dz^F = 0. It increases 
with temperature slower than in the PS model, due to 
the additional (overtwist) energy penalty of denatura- 
tion. The loop-size distribution essentially decays expo- 
nentially with a power-law correction as in the PS model. 

In the high-temperature phase where the loop fugac- 
ity is fixed at its upper bound zi = 1/s, the above pic- 
ture is no longer valid. Setting now dF{zi,mi,mo) = 
yields a unique solution with TOq ^ (the extremal value 
Too = does not yield a minimum.) Therefore, a fi- 
nite fraction of the base pairs is located in a macroscopic 
loop. Setting dF/dnio = yields zi = 1/s, as expected 
in the high-temperature phase. Accordingly, the proba- 
bility distribution function for the microscopic loop sizes 
p{l) ^ is now scale invariant. This picture is in con- 
trast with the PS model, since no microscopic loop sur- 
vives for T > Tc- An analytical expression for the mass 
fraction in microloops and the macroloop can be obtained 
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FIG. 1: The fraction of the denatured DNA (soUd) and the contribution from the microscopic loops (dashed) for the supercoiling 
(left) and overtwisting (right) scenarios, as a function of temperature. The parameters s = 5, = 3, c = 3.5, ^4 = 0.1 are the 
same for both figures, while the stiffness parameters for overtwisting (k = 1.0) and supercoiling (see Ref. [3) are comparable. 



from the remaining two minimization conditions wrt nii 
and zf. 

dF 



= 



dmi 



ozi 



- (1 + ACc) = 
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niQ + R 



mi 



l + R ' 



(7) 
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where R \ 



i+AQ ^ temperature-independent constant. 
The transition is continuous, since setting toq = (T — > 
r+) recovers the critical value of m; in Eq.(IS]) found from 
the low-temperature limit T T~ . Solving for toq in 
Eqs.([7M8]), then substituting Inw = —eb/kBT and using 
Eq.(|4l) one finds that the macroscopic loop fraction has 
the exact form 

mo = 1 ^ =, (9) 

where B = ^^(1 + i?)^^. mo grows linearly with tem- 
perature in the vicinity of Tc and approaches unity only 
in the limit T — >■ oo. The total fraction of microscopic 
loops decreases with temperature while their distibution 
remains a power law. In Fig IT] we present a comparison of 
the "overtwist" and the "supercoil" dominated scenarios 
in terms of macroscopic and microscopic loop fractions, 
for a generic set of parameters at which both display a 
second-order melting transition. Note that, the discon- 
tinuity in dmi/dT aX T = (inset) in the supercoiling 
model is hardly visible. However, the cusp in the mi- 
croscopic loop fraction is clear and exists even for c < 3 
where the discontinuity at shifts to higher derivatives 
of the free energy (see Ref . fio'] for more on PS model with 
supercoiling). The cusp in mi is more prominant in the 
"overtwist" picture with a comparable set of parameters. 

The presence of the macroloop above Tc does not de- 
pend on the precise value of c (as long as there's a melting 



transition), while its size relative to the total amount of 
denatured DNA, mo/mi, does (Figl2]). For c < 2, the 
macroscopic loop vanishes together with the transition 
itself. Otherwise, the relative size of the macroloop grows 
faster with temperature as c gets larger. This trend as a 
function of c converges to a limiting curve which is indis- 
tinguishable from that obtained for c = 3.5 and shown in 

Figia 

Generalization of our results to cr = {Lg — i;)/i 7^ is 
straightforward after substituting f3K{mi + o')^/(l — mi) 
for the twist penalty in Eq.(|3]). The melting picture for 
nonzero a remains qualitatively unaltered, except when 




(T-T^)/T^ 



FIG. 2: The weight of the macroloop among the total amount 
of denatured bases is shown as a function of temperature, 
for various values of c. The uppermost curve for c — 3.5 is 
indistinguishable from those obtained for higher values of c. 
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a < —m'[, where the critical point itself disappears and 
a macroscopic loop exists at all temperatures. 

Finally, noting that the limit k — recovers the first- 
order transition, it is interesting to see how the sharp 
denaturation in the PS model smoothens out with the 
introduction of twist penalty. This crossover is shown 
in FigIS] where the total loop fraction given in Fig[l] is 
extended into the K-dimension. The region between the 
two sheets that join at the critical line t = is the macro- 
scopic loop fraction which smoothly approaches unity as 
At ^ 0. 




FIG. 3: (Color) The loop fraction m of the DNA as a function 
of the reduced temperature t = {T — Tc)/Tc and the overtwist 
penalty k. The red surface shows the contribution of the 
microscopic loops to the total denatured DNA. The size of the 
macroscopic loop is the vertical distance between the upper 
(gray onhne) and the lower (red online) surfaces for t > 0. The 
back panel of the bounding box corresponds to the PS model. 
Note that, the reduced temperature t does not correspond to 
a single temperature value, since Tc is a function of k. 



To summarize, we reconsidered the melting thermody- 
namics of a DNA chain with fixed linking number. We 
assumed that the denatured loops appear by transferring 
LK to duplex regions through overtwisting. We showed 
that, despite the different melting scenarios observed in 
supercoiling and overtwisting pictures, a feature common 
to both is the appearance of a macroscopic loop which 
grows monotonously with temperature. While the total 
fraction of denatured pairs increases with temperature, 
the DNA mass in the microscopic loops decreases above 
Tc- This condensation phenomenon is analogous to BEC, 
except it takes place at high temperature. Whether it is 
dynamically accessible is an interesting question, since 
merging microscopic loops towards a macroloop entails 
diffusing denaturation bubbles across torsionally strained 
duplex regions. Investigations in this direction, as well as 
towards a joint theoretical framework that incorporates 
both overtwist and supercoiling are in progress. 
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